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ABSTRACT 


The  inviscid  flow  field  about  a  right  circular  cone  in  unsteady  planar  flight  is 

analyzed  by  a  perturbation  technique  which  is  an  extension  of  Stone's  treatment  of 

the  cone  at  small  yaw.  A  solution  is  found  in  the  form  of  infinite  series  in  the  time 

rates  of  change  of  the  pitch  rate  and  angle  of  attack.  The  linear  stability  derivatives, 

Cm  and  Cm>  as  well  as  "higher  order"  stability  derivatives  such  as  Cm.  andC 
q  a  q  a 

are  presented  for  a  wide  range  of  cone  angles  and  Mach  numbers. 


The  stability  derivatives,  Cm  and  Cm  ,  as  obtained  from  this  solution  are 

q  Of 

compared  to  results  obtained  from  second  order  potential  theory,  Newtonian  impact 
theory,  and  an  unsteady  flow  theory  due  to  Zartarian,  Hsu,  and  Ashley.  Both  the 
potential  theory  and  the  impact  theory  predict  that  C  rapidly  approaches  zero  at 


a 


high  Mach  numbers  while  the  present  theory  indicates  that  Cm  approaches  a  value 
which  is  on  the  order  of  10  percent  to  20  percent  of  C 

mq 


Numerical  results  obtained  from  the  present  theory  are  also  compared  to  ground- 
test  data  on  Cm  +  Cm. .  The  agreement  is  found  to  be  generally  good,  although  the 

q  a 

data  in  some  instances  indicate  a  pronounced  Reynolds  number  effect. 


The  numerical  results  for  the  "higher  order"  coefficients  are  used  to  predict  the 

effect  of  reduced  frequency  on  the  parameter  C  +  C  as  obtained  by  the  forced 

mm. 

q  Ot 

iii 


J 


oscillation  testing  techniques.  It  is  found  that  the  predicted  effect  is  very  small  over 
the  range  of  reduced  frequencies  likely  to  be  encountered. 
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1.  INTRODUCTION 


As  a  re-entry  vehicle  penetrates  the  atmosphere,  its  pitch  rate  (q)  and  angle  of 
attack  (a)  vary  in  an  oscillatory  manner  and  cause  the  flow  field  about  the  vehicle  to  be 
in  an  unsteady  state.  For  most  purposes,  q  and  the  rates  of  change  of  q,  and  Or  are 
small  enough  that  the  differences  between  the  unsteady  field  and  a  quasi-steady  field 
may  be  considered  negligible.  An  exception  to  this  rule  occurs  in  the  evaluation  of 
aerodynamic  forces  and  moments.  The  small  contributions  of  unsteady  effects  to  the 
normal  force  and  pitching  moment  appear  as  damping  coefficients  in  the  equations  gov¬ 
erning  the  rigid  body  motion  of  the  vehicle  and  play  an  important  role  in  determining  the 
loads  the  vehicle  must  withstand. 

In  this  work,  the  unsteady  flow  field  about  a  vehicle  (specifically,  a  right  circular 
cone)  is  examined  under  the  assumption  that  the  motion  of  the  vehicle  is  given.  The 
vehicle  motio..  appears  explicitly  in  the  equations  and  boundary  conditions  for  deter¬ 
mining  the  flow  field  in  a  body-fixed  coordinate  system  (Section  4. 1).  In  order  to  make 
the  boundary  value  problem  tractable,  perturbation  parameters  which  are  measures  of 
the  departure  of  »he  flow  field  from  a  steady-state,  axis-symmetric  field  are  introduced 
(Equation  5. 14)  anJ  used  to  effect  a  linearization  c.r  the  problem  (Section  5.2).  The 
first-order  effects  of  these  parameters  are  considered,  and  a  formal  solution  in  the 
form  of  an  infinite  series  which  gives  the  flow  field  in  terms  of  the  vehicle  motion,  is 
assumed  (Equation  5.26).  Thus  the  mathematical  problem  is  reduced  to  one  of  solving 
an  infinite  sequence  of  linear,  ordinary  differential  equations  (Equations  5.27)  with 
variable  coefficients  which  are  determined  by  the  steady  field  at  zero  yaw.  These 
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equations  are  solved  by  numerical  methods  (Section  6).  The  numerical  results,  when 
properly  combined  (Appendix  I) ,  give  the  unsteady  flow  field  to  the  first  order  of  the 
perturbation  parameters.  The  unsteady  flow  effects  on  the  normal  force  and  pitching 
moment  coefficients  are  obtained  by  integration  of  the  surface  pressures  (Appendix  III) . 

Results  obtained  by  the  method  presented  here  are  compared  to  theoretical  results 
from  other  methods  of  varying  degrees  of  approximation,  and  are  also  compared  to 
experimental  results  (Section  7).  From  the  cases  examined,  it  does  not  appear  that 
there  is  an  orderly  pattern  of  agreement  between  the  results  obtained  here  and  the  po¬ 
tential  theory  of  Tobak  and  Wehrend  (Reference  13),  the  shock  expansion  theory  of 
Z  art  ail  an ,  Hsu  and  Ashley  (Reference  14),  or  Newtonian  impact  theory.  Comparisons 
with  experimental  data  do  show  fair  agreement. 


Manuscript  released  by  the  author  September  1964  for  publication  as  an  FDL  Technical 
Documentary  Report. 
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2.  ASSUMPTIONS  AND  RESTRICTIONS 

The  flow  field  is  determined  by  the  body  geometry,  the  nature  of  the  gas,  and  the 
flight  conditions.  The  assumptions  and  restrictions  pertaning  to  each  are  tabulated 
below: 

1 .  Body  Geometry  -  The  body  is  assumed  to  be  a  right  circular  oone. 

2.  Nature  of  the  Gas  -  The  gas  is  assumed  to  be  inviscid,  non-conducting  and  at 
chemical  and  thermodynamic  equilibrium.  It  is  represented  analytically  as  a 
y*  gas  as  described  in  Reference  1 .  This  includes  an  ideal  gas  with  constant 
specific  heats  as  a  special  case. 

3.  Flight  Conditions  -  The  vehicle  trajectory  is  assumed  to  be  planar  with  three 
degrees  of  freedom  -  i.e. ,  two -degrees -of -freedom  in  translation  in  the  plane 
of  the  trajectory  and  one-degree -of -freedom  in  rotation  normal  to  the  place  of 
the  trajectory.  ft  is  further  assumed  that  the  speed  of  the  vehicle  is  constant 
and  sufficiently  high  that  the  flow  about  the  vehicle  is  supersonic  with  respect 
to  the  vehicle.  This  reduces  the  number  of  degrees  -of-freedom  to  two -angle 
of  attack  and  pitch.  Only  first  order  effects  of  these  on  the  flow  field  are 
considered. 

The  variation  of  ambient  pressure  and  density  along  the  trajectory  is  neglected  in 
the  analysis. 
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3.  DESCRIPTION  OF  THE  PERTURBATION  SCHEME 


The  flow  field  about  a  body  to  flight  is  determined  by  the  solution  to  the  non-linear 
boundary  value  problem  given  in  Section  4.1.  This  boundary  value  problem  is  stated  in 
a  coordinate  system  fixed  in  the  body  and  consequently  the  motions  of  the  body  $  and 


eg 

Q  ,  appear  as  "driving  functions"  in  the  problem.  For  the  simple  planar  trajectory 
considered  here  the  motions  of  the  body  are  given  in  terms  of  two  functions  of  time. 


Or  ft)  and  q  (t),  and  the  constant  speed  UB  =  |  | ,  by  Equation  (4.13).  The  flow  field 

variables  are  functionals  of  the  functions  o  (t)  and  q  (t)  in  that  they  depend  on  all  the 
values  taken  on  by  o(t)  and  q  (t)  in  the  interval  from  the  initial  time,  to  the  current 
time,  t,  and  are  ordinary  functions  of  position  as  given  by  the  three  coordinates. 


(R,  u\  <p) . 


fa  the  perturbation  scheme  utilized  in  this  work ,  or  (t)  and  q  ^.L  are  taken  to  be 

U» 

small  quantities  on  the  order  of  and  C2»  respectively,  which  then  become  the  pertur¬ 
bation  parameters  with  the  substitutions: 


o  <t)  *  a  (t) 

U. 

q  (t)  =  C2  q  (t)  ~ 

the  pressure,  for  example,  will  be  a  functional  of  the  functions  or  (t)  and  q  (t)  and  an 
ordinary  function  of  the  parameters  r  j  and  f  ^ .  It  is  assumed  that  p  can  be  expanded  in 
Taylor's  series  in  the  parameters,  (j  and  C2.  to  give  a  series  of  the  form 

P  =  Po  *  C1  pl  *  *2  p2  *  •  •  •  (H.O.T.) 
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The  coefficient,  pQ,  is  the  pressure  field  produced  by  the  body  in  steady  flight  at  zero 
angle  of  attack  and  can  be  found  by  established  methods.  The  coefficients,  p^  and  Pg, 
give  the  first  order  effects  of  angle  of  attack  and  pitching  rate,  respectively,  and  are 
to  be  determined  by  solution  of  Equations  (4.16)  and  (4.17).  They  are  functionals  of  the 
functions  a  (t)  and  q  (t)  and  ordinary  functions  of  the  spatial  coordinates.  In  Section  4.2 
it  is  shown  that  they  can  be  represented  formally  by  series  of  the  type: 


(and  a  similar  series  for  p2  as  given  in  Appendix  1)  where  the  pl  n  are  functions  only  of 
the  ray  angle,  co.  This  solution  holds  after  "starting  transients"  have  died  out.  The 
coefficient,  Pj.q.  gives  the  effect  of  small  yaw  in  the  steady  state.  The  coefficients, 
P|  ,j ,  pj  2 •  etc . ,  give  the  effects  of  time  varying  angle  of  attack. 

It  is  the  coefficients,  ps  (and  corresponding  quantities  for  the  other  flow  field 
variables) ,  which  are  found  as  a  result  of  this  analysis.  The  method  of  finding  them  is 
numerical  and  is  described  in  Sections  4  and  5.  The  pressure  coefficients  yield  cor¬ 
responding  forces  and  moment  coefficients  (Appendix  III)  which  are  static  and  dynamic 


stability  derivatives. 


4.  DERIVATION  OF  THE  PERTURBATION  EQUATIONS 


4.1  STATEMENT  OF  THE  BOUNDARY  VALUE  PROBLEM  IN  GENERAL  FORM 


The  inviscid  flow  field  boundary  value  problem  can  be  stated  to  an  observer  in  a 
body -fixed  coordinate  system  (x,y,z),  Figure  1,  by  a  transformation  of  coordinates 
from  an  inertial  system.  Lamb.  (Reference  3).  gives  the  appropriate  transformed 
continuity  and  Euler  Equations: 

the  continuity  Equation : 

dP  +  v.  pV  -  0; 

3t  ’  (4.1) 


the  Euler  Equations: 


dV  -*  nr  -*  ?p 

-2— +  (V.  0)  V  ♦  2  (Ox  V)  +  — 

at  o 


diT 


eg 


■*  -±  dfi  >  _  j. 

-  n  x  V  ♦  r  x-j£-*  ( <">x  r)  x  n. 


dt 


eg 


(4.2) 


The  form  of  the  continuity  equation  is  unaltered  by  the  transformation.  The  form 
of  the  Euler  equations  is  altered  in  that  the  acceleration  of  a  fixed  point  in  the  moving 
frame  of  reference  appears  as  a  body  force  (per  unit  mass)  involving  the  vector  v  elocity  of 
the  center  of  gravity,  V  (t),  and  the  vector  rotation  of  the  body,  0  (t),  on  the  right  hand 
side.  The  vector  rotation  causes  a  given  fluid  particle  to  experience  a  Coriolis  accel¬ 
eration,  2  (0  x  V),  which  is  present  on  the  left  hand  side.  The  derivatives,  ^Vcg 

dt 

and  appearing  in  equation  (4.2)  are  those  observed  in  the  body-fixed  coordinates. 

^  dV 

For  example,  is 

dt 
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SHOCK,  F  =  0 

'  R 


Figure  1.  In  vise  id  flow  Field  Boundary  Value  Problem 


^cg  -  d  <*  *  VCg) 
dt  '*  dt 


d  (y  •  *cg> 
dt 


+ 


d  <*  •  %£ 

dt 


The  continuity  and  Euler  equations  must  be  complemented  by  the  energy  equation 
for  the  adiabatic  flow: 


~  *  V  .  VS  =  0  (4.3) 

dt 


and  an  equation  of  state  for  the  gas,  here  taken  to  be: 


The  Auction  Indicated  in  Equation  (4.4)  depends  on  whether  the  gas  is  an  ideal  gas, 
a  gas  at  chemical  equilibrium,  or  a  gas  in  some  "frozen"  composition.  The  development 
of  the  equations  can,  however,  proceed  without  specification  of  the  precise  form  of  the 

equation  of  state. 

Equations  (4.1)  through  (4.4)  are  a  complete  set  of  flow  equations  for  the  determina¬ 
tion  of  the  pressure,  p,  the  density,  p,  the  entropy,  S,  and  the  three  components  of  the 
fluid  velocity  (measured  in  the  moving  frame  of  reference),  V,  for  prescribed  motions 
of  the  body,  V  and  0.  In  order  to  solve  them,  initial  conditions  at  some  instant  of  time 
and  botsidary  conditions  at  the  body  and  shock  must  be  given. 

As  initial  conditions  it  is  assumed  that  it  t  =  0,  the  field  is  the  steady-state,  axi- 
symmetric  field  produced  by  a  uniform  forward  translation  at  speed  U„  of  the  body 
along  its  axis  of  symmetry. 

At  the  body  surface,  the  flow  must  be  tangent  to  the  body  surface  ano  the  boundary 
condition  is: 

V  •  nB  »  0  (4.5) 

On  the  shock  surface : 

Fs  (x,;  *,t)  =  0,  (4.6) 

(which  must  be  found  as  s  part  of  the  solution)  the  shock  equations  give  the  flow  variables 
as  functions  of  the  relative  velocity  of  the  shock  and  the  free  stream  flow,  and  the  in¬ 
stantaneous  unit  inward  normal  to  the  shock. 

8 


i 


J  , 


<KHi 


"n  =*7  Fs 

8  |*Fj 

*4 

The  sign  in  equation  (4.7)  is  to  be  chosen  so  that  tc  .  n8  is  negative. 


The  shock  equations  can  be  reduced  to  the  following  set  of  three  algebraic 
equations . 


(4.7) 


ps  Vnn  ~  p®  vn« 

p  V2  +p=pV2+p  ^ 

g  ns  rg  *  u®  od 

\  2  y2 

ns  n06 

h  (p  p  )  +  - =  h  (p  0  )  ♦  - 

s,  8  2  ®  «  2 


These  equations  express,  respectively,  the  conservation  of  mass,  momentum 
normal  to  the  shock,  and  energy  in  a  form  which  utilizes  the  result  that  the  velocity 
component  tangent  to  the  shock  is  unchanged  in  crossing  the  shock.  In  these  equations 
Vn0D  is  the  component  normal  to  the  shock,  of  the  relative  velocity  between  the  shock 
and  the  flow  on  the  upstream  side  of  the  shock,  and  Vn8  is  the  corresponding  compo¬ 
nent  on  the  downstream  side  of  the  shock. 


The  flow  velocity  on  the  upstream  side  of  the  shock  is  given  by: 


^  -  V  *  {?**) 

•  eg 


(4.9) 


and  the  component  of  shock  velocity  along  its  normal  is  given  by  (reference  3,  page  7): 

d  F 
_ s_ 

it 
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Therefore,  VnJ,  is  given  by: 


Vnl  *  V« 


n8  + 


nQ  +  (r  x  0)  *  nR  + 


Equations  (4.8)  are  presumed  to  be  solved  in  the  form*: 

\ 


P8  =  P8  (Vn«.  P«.  P«) 

Pg  ~  Pg  (Yj) 00  ’  P®  >  P») 

Vn8  "  vn»  =  *  Vn  =  *Vn  <vn- P ..  P.)  , 


(4.10) 


(4.11) 


The  first  two  of  these  equations  give  the  pressure  and  density  downstream  of  the 

shock  explicitly  in  terms  of  the  function  giving  the  shock,  Fg,  and  the  motions  of  the 
«■ •  — ♦ 

body,  VCg  and  Oby  use  of  equation  (4.10).  The  third  equation  gives  the  three  velocity 
components  of  the  flow  on  the  downstream  side  of  the  shock  by  further  use  of  the  equation: 


Vs  = 


■V  + 

eg 


(r  x  0)  +  P  Vn  n8 


(4.12) 


4.2  REDUCTION  OF  THE  BOUNDARY  VALUE  PROBLEM  TO  PERTURBATION  FORM 


The  equations  and  conditions  of  Section  4.1  define  the  boundary  value  problem  for 

*"4 

determination  of  the  inviscid  flow  field  about  a  body  whose  motions,  V__  and  0  ,  are 

♦There  are  numerous  procedures  in  the  literature  for  obtaining  such  a  solution  numer¬ 
ically  .  Reference  4  presents  the  procedure  used  in  the  subsequent  numerical  work. 
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known.  In  this  Section,  the  boundary  value  problem  is  reduced  to  a  perturbation  form 

^  I 

applicable  to  right  circular  cones  which  move  at  constant  speed,  VCg  =  const.  = 
in  such  a  manner  that  the  z-axis  is  always  parallel  to  a  fixed  reference  direction.** 


For  the  type  of  motion  just  described ,  the  vector  velocity  and  rotation  can  be 
written  in  terms  of  two  functions,  a  (t)  and  q  <t): 


vCg =  u®("x  008  a+  y  8ina> 


n-  qz 


(4.13) 


To  put  the  problem  in  perturbation  form  it  is  assumed  that  a  «  1  and  ^ —  « 1, 

CO 

and  the  parameters  Cj  and  e2  which  measure  the  magnitudes  of  a  and  ,  are  intro- 
duced  by 


o  =  Of 


qL  - 

811(1  u  =  ^2  q 
00 


(4.14) 


In  Equation  (4.14),  a  (t)  and  q  (t)  are  assumed  to  be  of  the  order  of  one  in  the  time 
interval  of  interest,  and  and  c2  (which  are  constants)  are  assumed  to  be  much  less 
than  one.  The  flow  field  variables  can  be  regarded  as  functionals  of  the  functions  a  (t) 
q  (t)  (in  the  sense  defined  in  Reference  3)  and  ordinary  functions  of  the  perturbation 
parameters  and 


♦♦Although  this  analysis  does  not  employ  stability  axes  the  force  and  moment  derivatives 
obtained  are  conventional. 
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In  the  perturbation  scheme,  it  is  assumed  that  each  flow  field  variable  can  be  ex¬ 
panded  in  Taylor's  Series  in  the  variables  c  ^  and  c  2»  For  example,  it  is  assumed  that 
the  pressure  can  be  written  in  the  form: 


P  =  Po  (x,y,z)  +  Ci  Pi  (x.y.z.t)  +  e 2  P2  (x,y,z,t)  +  . .  .  .  H.O.T..  (4.15) 

These  expanded  forms  are  then  substituted  into  the  equations  defining  the  boundary 
value  problem  and  the  resulting  set  of  perturbed  equations  is  separated  into  a  number 
of  sets  of  equations  for  the  perturbations  by  the  usual  process  of  equating  forms  in  like 
powers  of  and  C2*  (Only  the  three  lowest  order  terms  are  considered.) 

The  boundary  value  problem  so  defined  for  the  lowest  order  coefficients,  P0»P0»S0 
and  VQ,  is,  of  course,  the  problem  of  determining  the  steady-state,  axisymmetric  field 
produced  by  the  body  as  it  translates  with  constant  speed,  Um ,  parallel  to  its  axis  of 
symmetry.  Under  the  assumptions  listed  in  Section  2,  this  field  will  also  be  conical 
and  can  be  solved  by  the  methods  of  Taylor  and  Maccoll  (Reference  5) ,  or  Romig  (Ref¬ 
erence  6).  In  the  subsequent  work  it  is  assumed  that  this  field  is  known*.  The  equations 
defining  the  steady  state  field  are  given  in  Appendix  n. 

Substitution  of  Equation  (4.15)  and  like  expressions  for  the  density,  entropy,  and 
fluid  velocity  into  Equations  (4,1)  through  (4.4)  to  obtain  the  two  sets  of  equations  gov¬ 
erning  the  coefficients  of  Cj  and  yields. 


♦The  numerical  procedure  used  to  obtain  it  is  given  in  Reference  7 . 
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i 


dp. 


d  t 


av 

d 


I  -  -  -  -  7  Pi  ?P0 

-  +  (V- *)  V  +  (V4.  v)  v  +  -  -2  p  =  F 

t  0  1  3  0  2  J  J ; 

P0  P0 


\ 


(4.16) 

J  =  l,2 


dS]  - 

- L+  V  •  VS  +  V  •  vS  =  0  ; 

at  0  i  J  o 


pj  =  co2  pj  +eoSj- 


where: 


-  U 


da 


dt 


1=1 


V 


U  U  Udq 

2  (V  x  z  )  q - +  —  qy  +  (rxz)  —  -  , 

L  L  L  dt 


J  “  2 


(4.17) 


The  boundary  condition  at  the  body  surface.  Equation  (4.5),  separates  into  the  two 
sets  of  conditions: 

Vj  *nB  =  0,  0  =  1,2)  (4.18) 

For  the  conical  geometry,  the  equation  giving  the  shock,  Equation  (4.6),  is  most 
conveniently  expressed  in  terms  of  the  spherical  coordinates,  (R,  to,  <p),  Figure  2,  in 
a  form  giving  the  ray  angle,  to,  as  an  explicit  function  of  the  radial  coordinate,  R, 
meridional  angle,  <p,  and  time,  t: 

Fs  Ha's  (R,  <p,  t) -to=  0  (4.19) 
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Figure  2.  The  Spherical  Coordinate  System  (R,  co,  <p) 

The  perturbation  form  for  the  shock  is  then  (to  the  three  lowest  order  terms) : 


cc  =  U's  -  c*fe0  +  ( l  a'sl  <R»  <P>  ^  +  €2  ws2  (R’ <p*  ^ 


(4.20) 


where  u  is  the  constant  shock  angle  from  the  steady -state,  axisymmetric ,  comical 
flow  solutions.  The  unit  inward  normal  to  the  shock  (Equation  4.7)  becomes: 

_  (  _  dw  .  cp  d  uu  j 

;  \  c.  )r  r_±  ♦  _  J  I 

8 


E 


dw  .  cp  d  w 

R  R  _!>  *  _ 


SR  sin  uu  dcp 

so 


(4.21) 


and  the  flow  velocity  on  the  upstream  side  of  the  perturbed  shock  is 
V  -  R  f  u  _  ♦  c,  u  .  +  e„  u_A  +  v  ♦  e  .v  ♦  e, 


■  "  |  U.o  *  £l”.l  *  C2  ".zj*  *{V.o  *  'lV.l  * 

®{cl%  *  *2  *.2} 


V-2  K 


(4.22) 
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where: 


u  =  U  cos  UU 
®  O  w  SO 


U  ,  =  4j(a 

®  1  »i 


sin  uu  cos  rp  +  uu 

so 


,  sinuu  J 
si  sO  f 


u  =  x  gin  uj  cos  rp  -  U  uu  .  sin  uu 

®  2  T  eg  so  ®  so 


v  =  -  U  sin  uu 

» O  w  SO 


(4.23) 


v  =  -  U  (a  COS  UU  CO8C0  +  uu  ,  cos  uu  ) 
®1  ®'  80  ^  si  SO  ' 


v»2  = 


_qUa 


(  x  COS  UU  -  Rj  COP  tp-  U  0)  „  COS  UU 
Cfc  SO  ®  s?  so 


w  =  a  U  sin  <n 
00  1  00 


w 


®  2 


«UW/  \ 

- IRcoSuu  -x  1  sin  cp 

L  V  80  Cg/ 


The  perturbation  form  for  the  normal  component  of  relative  velocity  between  the 
upstream  flow  and  the  shock  is  obtained  from  Equation  (4.10)  as: 


V  =  -  V  + 
n®  cd  o 


mt 

£ 


c.  v , 

1  nj 


J  =  1 


where  t 


(4.24) 


duu. 


duu 


V  ..  R_!i  *  RZfi  -  V 

dR  *t 


Substitution  of  Equation  (4.24)  into  Equations  (4.11)  and  (4.12)  gives  the  pertur¬ 
bation  forms  of  the  pressure,  density,  and  velocity  components  on  the  downstream  side 
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of  the  perturbed  shock.  These  resulting  expressions  must  be  equated  to  expressions 
for  the  corresponding  variables  in  the  field  evaluated  at  the  perturbed  shock  to  obtain 
the  proper  boundary  conditions  for  the  perturbation  equations.  The  resulting  condi¬ 
tions  are:  (j  3  1, 2  in  each  equation) 


PJ  ^R*  u'so’  cp*  ** 


dp 


d  uu 


3P«, 


UU  .  + 

Sj 


av. 


UU  =  U) 


n® 


so 


nj 


V  =  -  v 
n®  oo  o 


1 


(R,  tv  ,  qp,  t) 
so 


duu 


- 


uu 


8j 


av, 


Ul~  IJU 


n«o 


so 


nj 


V  =  -  v 
n®  ®o 


V  (R,  uu  cp,t)*(u  R  +  v  ,  uu  +  w  cp)-( 
j  so,  ®j  ®  j  ®j 


'd  u 


uu 


R 


L)  =  uu 


dv  ^ 
°>w 


du) 


so 


uu 


8j 


UU=  uu 


so 


auu  i  i  duu 

+  A  V  (  -  v  pp)(R  SJR  1  ®J5P 

n  ®0,  a  ,  I  -  f 


SR 


sin  ai  d  cp 
so 


dA  V 


dV 


n  - 
V  .  x 
nj 


n® 


V  =  -  v 
n®  x  o 


(4. 25) 


In  Equations  (4.25),  the  partial  derivatives  with  respect  to  u:  arise  as  a  result  of 
evaluating  the  variables  in  the  field  at  the  perturbed  shock.  The  partial  derivatives  with 
respect  to  Vnj  are  to  be  obtained  by  differentiation  of  the  expressions  indicated  in 
Equation  (4.11);  expressions  for  them  are  given  in  Reference  9. 
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To  complete  specification  of  the  boundary  value  problem  for  the  perturbation  vari¬ 
ables,  it  is  necessary  to  give  the  initial  conditions  at  t  =  0,  Corresponding  to  the 
assumption  stated  in  Section  4,1,  this  condition  is  that  the  perturbation  variables  vanish 
at  t  =  0. 


Equations  (4.16)  and  (4.17)  together  with  the  boundary  conditions  at  the  body  sur¬ 
face,  Equation  (4.18),  the  boundary  conditions  at  the  shock,  Equation  (4.25),  and  the 
initial  conditions  define  two  separate,  linear  boundary  value  problems  for  the  perturba¬ 
tion  variables.  In  these  problems,  the  "motions"  of  the  body,  a  (t)  and  q  (t),  appear  in 
inhomogeneous  terms  and  play  the  role  of  'driving  functions. "  Because  of  the  linearity 
of  the  problems,  the  solutions  for  arbitrarily  prescribed  functions,  tv  and  q,  can  be 
divided  into  parts  which  can  be  called  the  particular  solution  and  the  complementary 
solution.  The  ^articular  solutions  are  defined  as  solutions  of  the  inhomogeneous  prob¬ 
lems  which  reduce  to  the  trivial  solution  (all  perturbation  variables  equal  zero)  when  a 
and  q  are  Identically  zero;  the  complementary  solutions  are  defined  as  solutions  of  the 
homogeneous  problems  (obtained  by  setting  a  and  q  equal  to  zero)  which  cause  the  com¬ 
plete  solution  to  satisfy  the  initial  conditions. 

On  physical  grounds,  it  is  known  that  the  complementary  solutions  must  'die  out" 
in  time,  and  their  decay  time  is  some  nominal  multiple  of  the  time  it  takes  the  body  to 
move  through  its  own  length,  L/U®  .  Subsequent  to  this  decay  time,  the  flo«  field  about 
the  body  is  given  by  the  particular  solutions  alone.  It  is,  therefore,  only  the  particular 
solutions  which  are  of  interest  here. 
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The  problem  of  determining  the  particular  solutions  of  the  boundary  value  prob¬ 
lems  for  arbitrarily  prescribed  functions,  a  (t)  and  q  (t),  can  be  approached  in  several 
ways  —  e.g. ,  by  use  of  the  Laplace  transform.  The  approach  taken  here  is  one  which 
reduces  each  problem  to  a  problem  of  solving  an  infinite  sequence  of  sets  of  ordinary 
differential  equations  having  the  ray  angle,  u>,  as  an  independent  variable.  The  substi¬ 
tutions  which  accomplish  this  reduction  are  tabulated  in  Appendix  I.  The  form  of  the 
pressure  perturbation,  plt  is  typical: 


<R»  u» .  «p .  t )  ■  \  '  p  b»)f  \  d  *  cos  cp 

^  •  W 

n  =  o 


(4. 26*) 


In  writing  these  series,  it  is  assumed  that  a  (t)  and  q~ (t)  are  analytic  functions.  Also, 
it  is  recognized  that  there  is  no  a-priori  guarantee  that  these  series  are  convergent  for 
any  given  functions,  S  (t)  and  q  <t).  The  factor  l/u”  has  bee  introduced  so  that  the 
perturbation  coefficients,  p.  have  the  units  of  pressure. 

The  two  infinite  sequences  of  sets  of  ordinary  differential  equations  for  the  varia¬ 
bles  Pj  q(jj),  etc. ,  which  arise  from  the  substitutions  listed  in  Appendix  I  can  be 
written  in  the  common  form,  (j  »  1,2),  (n  »  0, 1,2,3.  ... ). 


j(|Do'',,o**oOi.n) 

*  Pv  wj,n  “^J.n 


stn  us 


(4. 27a) 


•The  "zeroth”  derivative, -  ,  is  equal  to  a,  by  convention. 

dt® 
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(4.27b) 


^J.n 

vo  - v°  V 


J,n  +  <n+*],2>  (uouJ.n+^)  ’^J.n 


dv 


IlB 


do: 


K  1  !ii»  dp«  , 

♦“o’j,n*vo“i,n*v),n -  * - - - 2  Z 

du;  Pc  dw  PQl  dw 


(n  +  6j  2)  “ovj.na4.' 


(4.27c) 


dw. 


du; 


+  u. 


wj»n 


v0  w. 


J.o 


cot  u; 


-it9 -  +  (D  +  6.  )|L,«(  .  n 

>.  sinu;  '  J.2  ^  J'n  J*n 


(4. 27d) 


v 


o 


Sla» 

du.' 


+  (“  +  6j,2)uoSj,n“-^,n 


(4.276) 


(4,  27 1) 


where: 


^i.n* 

-(1 

A.»' 

-(1 

-  6o,o>  U-  “l.n-1  *  6l.o  4J.  1  u2-  •*"“*  Vo  6).2  u-  <2  vo  * 

A.- 

-(1 

-  «o.n)  «-  (*i.l  co»  «*  V*>  * 

6o,n 

6 2  U.  (-2  Uq  ♦  U.  cos  uj) 

-(1 

-  V«>  U-  wi.n-l  *  4l,n  U.2  <»J.l  *  !),2  co*  “)  * 

6o,n 

6j  2  j2  U.  (uQ  co.  ic  -  v0  sinu;)  -  U»2) 

x^j.n  = 

-(1 

-Vo>u*s).o-l 

The  condition  of  irrotationality  of  the  field  at  aero  yaw  (Appendix  ID  baa  been  used 
to  aimplify  equation  (4. 27b). 


19 


The  boundary  conditions  for  these  equations  can  also  be  written  in  a  common 
form.  At  the  body,  u  •  a'g,  the  condition  of  tangency  of  the  flow  becomes: 

vjt n  •  0  at  <c  •  u:B  (4. 28) 

At  the  shock  a*  •  u;  ,  the  perturbation  variables,  p,  „  etc. ,  are  related  to  the 
so  j,n 

quantities,  u.  which  describes  the  perturbation  m  the  shock  shape  by: 


UJ, 

J.n 


(4.29a) 


dpo 

*P2 

| 

pj.n* 

3u; 

+  (n  +  1  +  6j  2)  U.  cosaso 

u:  ■  u;  ~vnl 

80 

V»1  =  *v-oj 

'ix. 

).n 


♦  \  » 


*02 


J.n  dV 


nl 


V  -  _v 
nl  xo 


(4.29b) 


(4.29c) 
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(4. 29d) 


Uj  n  =  (n+l  +  6J2)  ^Vn(-v.0|  P->  P»)  W},n  -60,n6Mu«  £in  ‘So 


K-  &h.. 


so 


SAV, 


-<"*1*6I,2>  U.  cosSo  W 


nl 


V  *  ■  -  v 
nx  ®  q 


l.n 


-  6. 


,  r  .  1  ,  **vn 

,nu-  l4J,l  »J.2j 


'nl  ~v»o> 


(4.  29e) 


W;  „ 

J.n 


AV  /-  v„  p  p  ) 
n  1  o,  ,  ® ' 

sin  u^o 


*j,n  +  6o,n  U«<6j,i  +  6j,2  cos  *So>  *4“29f> 


where: 


X.  =  6 
J.n  o.n 


L:  (6.  cos  i*.  ♦  6  )  +  (1-6  )  U 

®  j,  1  so  j,  2  o.n'  * 


i.n-l 


Equations  (4.27a)  through  (4.27f)  are  ordinary  differential  equations  with  variable 
coefficients  which  are  determined  by  the  solution  ior  the  axially  symmetric  conical 
field.  u0,  v0>  p0  and  pQ.  For  j  *  1  and  n  •  0,  the  Equations  are  equivalent  to  those  de¬ 
rived  by  Stone  (Reference  9)  for  the  problem  of  cones  at  small  vaw.  For  eithei  value 
of  j  and  any  value  of  n  (except  n  -*  0),  the  equations  and  boundary  conditions  are  com¬ 
plete  provided  that  the  solution  for  the  same  value  of  )  and  a  value  of  n  which  is  one 
less,  has  been  previously  solved.  For  n  «  0,  the  equations  and  boundary  conditions 
are  complete  without  a  knowledge  of  the  solution  for  any  other  value  of  j  and  n.  For 
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any  values  of  j  and  n,  the  problem  is  a  two-point  boundary  value  problem  in  which  the 
boundary  conditions  contain  an  unknown  constant 

The  equations  have  a  singular  behavior  at  the  body  surface  where  v0  =  0  since  v0 
multiplies  the  derivatives  of  the  perturbation  variables  in  equations  (4.27b)  through 
(4. 27e).  The  method  of  treating  this  singularity  is  given  in  Section  5. 


5. " 


K » 
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5.  NUMERICAL  SOLUTION  OF  THE  PROBLEM 


The  numerical  method  of  solution  of  the  problem  is  based  on  manipulated  forms 
of  Equations  (4.27a)  through  4.27f).  The  objective  of  these  manipulations  is  to  reduce 
as  much  as  possible  of  the  numerical  work  to  performing  quadratures  and  also  to  pro¬ 
vide  a  convenient  method  of  handling  the  singularity  at  the  body  surface. 

The  manipulated  forms  of  the  equations  are  obtained  as  follows: 


Dividing  through  equation  (4.27e)  by  vq  and  solving  it  by  use  of  an  integrating 
factor  gives 


U) 


so  so 


where 


and 


“  ■  ”  +  {2,J 


tli  =J-  P  v  sin  uo 

Of  0  0 


has  been  introduced 
v 


o 


u  = 
0 


by  use  of  the  equation 


d  \b 

o 

d 

co 


from  Appendix  n. 


23 
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Multiplying  equation  (4.27b)  by  u^,  equation  (4.27c)  by  v  ^  and  equation  (4.27e)  by 

T  ,  adding  the  results  and  combining  terms  gives 
o 


dH. 

v  — — +  mu  H.  =  u  .  +  v  ^  +  T  */ 

o  d  o  ),b  o  j,n  o  j,n  o  j,n 

LC 


(5.2) 


where  H.  is  a  quantity  related  to  perturbations  in  the  stagnation  enthalpy 

J.n 

Pi  n 

H  =?  +  u  u.  +  v  v.  +  T  S. 

j,n  p  o  j,n  o  j,n  o  j.n 


(5.3) 


The  condition  of  irrotationality, 


du 

c 

Vo  "  d 


has  been  used  In  deriving  equation  (5.2).  Also  a  quantity, 

dT  p  dp  p.  dp 
s  _ o  _  j,n  _o  ^  J,n  o 

j,n  d  "  2  dw  2  dw 

p  P 

o  Mo 

which  appears  in  the  derivation  is  shown  to  be  zero  by  use  of  equation  (4.27f),  the 
relation , 


^o 

SCO 


1 


2  dU!  ’ 


the  identity, 
3  T 


3  p/s  p2  c  2 


and  the  fact  that 


S  =  const, 
o 
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The  form  of  equation  (5.2)  is  the  same  as  equation  (4.27e)  and  the  solution  is 


X 


H 


H 


j.n] 


j,n  \  fm 
o 


Y  m  +  Y 
o 


u  +v  4.  +  T  J 

i .  n  ft  In  ft  1 


to  =  U) 


SO 


m  f  o  ^j.n  vo  ),n 

0  J  .  m 

vi i 
o  o 

UU 

so 


O  j,n 


duu 


(5.4) 


Solving  eq.  (5.4)  for  p.  ,  substituting  into  equation  (4.27b)  and  re-arranging  the 

3  »n 


result  gives 


j.n 


d  u 


.  ^ .  -  m  H.  +  m  T  S, 

J'n  .  (m+nv  -  J>n  3>n  £  J»n 

duu  '  j,n  v 


(5.5) 


Substituting  for  p,  from  equation  (5.4)  and  v.  from  equation  (5.5)  (in  terms 
J.n  J.n 

of  J  )  in  equation  (4.27d)  and  re-arranging  again  gives  a  result  of  the  same  form  as 

J.n 

equations  (4.27e)  and  (5.2).  The  solution  is 


u 

W,  =  w-  +  , - — 

j, n  j.n  (m  +  1)  sin  x 


Uj,  n’ 
w,  sin  x  +  — — - 
j,  n  m+1 


■if 


m+1 


m+1 


sin  x 


x  =  x 


so 


.  ■*!  f  ( ii  _  t  S  V°J)-"  ] 

f  -ii  A  j."  0  M  (m+1)  l 

~  /  ^  J,  n  sl"  w  +  ' - ^ - 6“ 

iJ  vo*o 
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(5.6) 


The  preceding  results  cm  be  used  to  eliminate  v  ,  p.  and  W  from 

j.n  j.n  J.n 

equation  (4.27a)  and  thus  reduce  it  to  an  equation  for  u.  .  Equation  (5.5)  gives  v. 

.  J.n  J.n 

S.*. 

in  terms  of  J,  and - and  equation  (5.4)  subsequently  gives  p.  in  terms  of  H 

J.n  d  ui  j.n  j.n 

u  ,  ^  Jj  n  and  S  .  Substituting  this  result  into  equation  (4.27f)  then  gives  p  in 
J.n  —  ■* —  j.n  j.n 

do. 
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f_  =  -  o  W.  -  (m+2)  sin  <ju 

2  j,n  o  j.n 


0  u 
o  o 


H 


j.n 


P  u  v  j 

ooo  j,n 
c2  (m+1) 


-  u 


(5.12) 


P  T  +  e 
0  0  o 


S 


j.n 


Equations  (5.1),  (5.4),  (5.5)  and  (5.6)  give  S.  ,  H.  ,  J.  and  W.  as  linear 

j.n  J.n  j.n  j.n 

functions  of  the  shock  perturbation  parameter,  to.  by  substitution  from  the  boundary 

j.n 

conditions.  Equations  (4.29).  That  is,  they  give  the  variables  in  the  forms: 

S  =  S(1)  +  Sf>  <i> 
j.n  j,  n  j.n  j.n 

H  =  +  H(2)  uj 

j.n  j.n  j.n  j.n 


J  =  J(1)  +  J(2)  U) 
j.n  j.n  j.n  J,n 


(5.13) 


and: 


W  =  W(1)  +  W<2>  uu 
j.n  j.n  j,  n  j,n 


I 


where  the  superscripted  variables  are  kncwn  functions  of  co  (or  are  known  constants  in 
some  instances)  which  can  be  determined  numerically  by  quadratures.  With  the  super¬ 
scripted  variables  in  Equation  (5.13)  known,  the  quantities,  r  and  T  ,  which  determine 

1  4» 

the  righthand  side  of  Equation  (5.7)  are  also  known  in  the  sense  that  the  superscripted 
quantities  in  the  expressions : 


r  -rW  t®., 
11  1  j.n 


r  =  r  ^  +  r  ^2)  uj 

2  2  2  j,n 


(5.14) 
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are  known  functions  of  co.  Thus,  Equation  (5.7)  Is  a  second-order  ordinary  differential 
equation  of  the  form: 


dJu 


d  w 


l.P  ,  1  dA 

2  A  d(u 


d  u 


IdL 


d  uu 


+  B  u 


J.n 


T  +  U), 

1  2  j,n 


(5.15) 


where  r  and  t9  are  known  functions  of  to,  and  co  Is  an  unknown  constant.  The 
12  j  ,n 

solution  of  this  Equation  Involves  two  additional  constants  of  Integration,  K  and  K 
and  can  be  written  in  the  form 


u  =  sW  +  .«  m  +  g,  K  +z  K 
J.n  J,  n  J.n  j.n  *J,n  1  j.n  2 


(5.16) 


where: 


z .  is  a  solution  of  the  homogeneous  equation , 

J.n 


.2 

d  z  ,  ..  dz 

-JJL  .  »  «  -Jil  .  B  «, 
.2  A  dm  du;  J,  n 

d  uu 


»  0 


Subject  to  the  conditions, 
d  z 


ZJ,  n  (ib'  =  °'  d. 


iii 


'2  u0  (i'  B>' 


(These  boundary  conditions  are  chosen  so  that  z.  „  =  v  ) ; 

1,0  o 

g  is  a  linearly  independent  solution  of  the  homogeneous  equation  (Reference  10, 

J.n 

PK  33), 


(5.17) 
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J 


and  u(1)  +  u<2)  u  Is  8  particular  solution  of  equation  (5.15),  (Reference  10, 
j,n  j,n  j.n 


page  30) 


u<l)  ♦  u<2)  B  -  *. 

j.n  j,n  j.n  j 


/Ag  (t ;  +  a.  „  r9)  da 
j»n  '1  i.«  2 


The  particular  solution  is  separated  into  the  equations 


■S  -v/  A'l.nTkd“-gi.«  /  AZi."Tt 


d  w  (k=l . 2) 


Putting 


Tk  = 


(m  +  i)  r_L 

TT“  \d  a 


5ln“ 


‘ r,") 


and  A  from  equation  (5.8)  into  the  expression  for  u.  ^  and  integrating  by  part 
d  T  (k) 

eliminate  - 1 - gives 

d  a 

■ft  -  } <sl 


«  — * 


# 
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The  equations,  as  written,  contain  several  limiting  forms  of  the  type  0.®  and 

Ik) 

«  -®as  the  body  surface  is  approached.  With  one  exception,  J;  '  ,  these  limiting 

J.n 

forms  have  finite  values  and  the  final  results  for  the  perturbation  variable  are  all 


finite  at  the  body  surface.  The  following  limiting  forms  are  encountered. 


1.  Equations  (5.1),  (5.4),  and  (5.6)  contain  terms  of  the  type 


uu 

I  -  ♦  m  /  — da  where  m  +  o  (for  m  =  0,  f  3  0  in  all  cases)  and 

o  /  ,m 

j  v  * 
uu  o  o 
so 


f  (uu)  is  non-singular  at  the  body  surface. 

ID 


0 


f  (a) 

m 

v  t 

o  o 


dx 


As  the  body  surface  is  approached, 


However,  use  of  L'Hopitals  rule  and 


f  (tig) 

the  equations  of  Appendix  n  gives  lim  I  =  - : — ,  .  These  terms  are 

mUo(V 

handled  numerically  by  use  of  quadrature  formula  having  the  same  singularity 
as  the  integral. 


2.  Equation  (5.5)  requires  division  by  zero  at  the  body  surface.  For  ra  -  0, 
the  numerator  of  the  ft  ■S.V.-n  Is  identically  ze:t >  so  the  limit  is  zero;  for 
other  values  ot  m,  the  limit  is  not  finite.  However,  this  infinity  does  not 
appear  in  the  subsequent  numerical  work  as  explained  in  items  4  and  5. 
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3 


.  Equation  (5.17)  for  g  approaches  the  form  0.®>at  the  body  surface  due  to 

J.n 

the  boundary  conditions  imposed  on  z.  .  The  limit  is  again  finite  and  is: 

J.n 


gJ.n*V  “  2uq^b)  A(uub)  ‘ 


The  first  derivative  of  g.  (which  is  needed  in  Equation  (5.18)  approaches 

J.n 

the  form  •  -  •  at  the  body  surface  and  the  finite  value  of  this  limit  is  ob¬ 


tained  as  follows: 


differentiating  g.  in  equation  (5.16)  gives 

J.n 


U' 


As  an  identity  we  have 


or,  using  the  equation  tor  : 

J.n 
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Substituting  this  result  into  the  expression  for  the  derivative  of  g. 

J.n 


gives 


u'  =  a,'  so 

GO 


which  is  an  equivalent  expression  whi.-h  does  not  approach  a  limiting 
form  at  the  body  surface. 


00 

4.  In  the  computation  of  u.  ,  Equation  (5.18),  integrals  of  the  form: 

J.n 


I  = 


it 

I 


J,  dx 
j.n 


so 


where  e  is  free  of  singularities,  are  required.  Although  J  is  infinite  at 

J.n 

the  body  surface,  these  integrals  are  finite.  Bv  substituting  fnrd  from 

J.n 

euuation  (4.27e)  in  equation  (5.2)  and  solving  the  result  for  £  ,  equation 

J.n 

(5,5)  can  be  written: 


d  H 


d  S. 


M  -  T  -i’.?  -  b 

d  i  o  di  j.n 


j.n 


dH 

I'sing  this  expression  for  .!  and  integrating  bv  parts  to  'emove 

dS  >•"  '  du' 

and 


J.n 
d  u: 


gives  the  foilowing  expression  for  the  integrals 
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I  -  - 


-7  T  S  - H  \ 

u  l  o  J,n  j,n j\ 


— -  (H,  -  T  S.  ) 

u  ),n  o  j, n 
o 


'll  -  !1) 


SO 


<JJ 


♦  f  S  )-H.  J-Puj-  Jill 

J  J.n  da  \uq  /  j,n  du)  \  of  uq 

u; 

SO 


duu 


There  is  no  problem  with  singularities  with  this  form  for  the  integrals. 


uu 

/ 


5.  At  the  body  surface,  the  term 


§duj 


„  <k)  so  .  (k) 
* ,  e  J. 

1  J,  n 


apperring  m 


A  m+1 

Equation  (5.19)  assumes  the  limiting  form  »  -  ®  (except  when  m  =  0).  How- 
ever,  when  tfc*  expression  for  arising  from  Equation  (5.11)  is  sub¬ 

stituted  algebraically  into  this  expression,  the  infinities  caicel  without  the 
application  of  L'Hopitals  rule  to  give  a  finite  result. 

The  numerical  solution  of  the  problem  proceeds  in  the  following  order  (the  axisym- 
metric  conical  flow  solution  is  assumed  to  have  been  pre -computed). 

(1)  (2) 

1.  Compute  S  '  and  S  '  from  the  formulae  which  arise  from  substituting 

1  j  V  1  f  v 

the  boundary  conditions  at  the  shock,  Equations  (4.29a),  (4.29b)  and  (4.29c) 
into  Equation  (5.1). 


2.  Compute  and  from  the  formulae  which  arise  from  substituting 

1  ^  v  1  J  V 

the  boundary  conditions  at  the  shock,  Equations  (4.29a)  through  (4.29e)  into 


Equation  (5.4). 
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3.  Compute  W^1  q  and  w|2q  from  the  formulae  which  arise  from  substituting  the 
boundary  conditions  at  the  shock,  Equations  (4. 29d)  and  (4. 29f),  into  Equa¬ 
tion  (5.6). 

d  Zf  o 

4.  Compute  Zj  q  and  - — as  described  following  Equation  (5. 16),  using 
finite  difference  methods  when  m  ^  0. 


dgl  0 

5.  Compute  q  and  —  J  from  Equation  (5. 17)  and  its  derivative. 


6.  Compute  1^,  lj,2\  1^,  and  1^  and  the  remaining  factors  present  in  Equa¬ 
tion  (5. 18). 


7.  Compute  v^1  q  and  v|2q  from  Equation  (5. 19a). 

8.  Compute  the  constants,  U)1>0,  Kj  and  K£  in  Equations  (5. 16)  and  (5. 19)  from 
the  boundary  conditions,  Equations  (4.28),  (4.29d),  and(4.29e). 


9.  Compute  uA  0  from  Equation  (5.16),  v10  from  Equation  (5.19)>  Sjq.  h1,0» 
and  Wj  ^  q  from  Equations  (5. 13),  w^  q  from  Equation  (5. 6),  p^  q  from  Equa¬ 
tion  (5.3)  and  Pj  q  from  Equation  (4.27f). 

10.  Repeat  the  preceding  sequence  for  n  =  1,  2,  3 .  .  .  Nmax  and  for  j  =  2, 
n  =  0,1,2  ..  .  Nmax  to  generate  as  many  perturbation  coefficients  as 
desired. 
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6.  SAMPLE  RESULTS  AND  COMPARISONS 


Figures  3  through  14*  show  distributions  of  the  perturbation  coefficients  defined  in 
the  preceding  Sections,  across  the  shock  layer  of  a  10-degree  cone  at  Mach  10  for  ideal 
gas  conditions.  These  results  show  a  tendency  for  the  magnitude  of  the  perturbation 
coefficients  to  decrease  with  increasing  n.  This  tendency  continues  through  n  -  10  (the 
highest  value  for  which  results  iiave  been  obtained)  and  supports  the  use  of  series  of 
the  type  given  in  Equation  (4.26)  in  obtaining  a  solution  to  the  unsteady  flow  problem. 

The  perturbation  coefficients  are  fairly  well  behaved.  Near  the  body,  some  -  such 

-1/2 

as  those  in  entropy,  Figures  11  and  12  -  have  derivatives  which  vary  like  (u;  -  u.  ) 

B 

due  to  the  vq  which  multiplies  the  derivatives  in  Equations  (4.27b)  through  (4.27e). 

This  behavior  does  not  cause  any  difficulty  in  the  numerical  work,  however. 


The  perturbation  coefficients  of  the  shock  angle  corresponding  to  the  results  shown 
in  Figures  3  through  14  are  given  in  the  table  below. 

TABLE  I 

PERTURBATION  COEFFICIENTS  FOR  THE  SHOCK  SHAPE 
ujb  =  10°,  M  x  -  ID,  >=1,4 


u.'  =  0. 215  rad 

so 


1,0  ' 

-  0. 0969 

Cl'  — 

2,0 

-  0. 176 

1,1 _ 

-  0.0742 

1 

-  0.0112 

1.2  ' 

0.  0319 

CL 

2  2- 

0.  0104 

1,3  ' 

-  0.  00737 

4-,  J 

-  0.  00250 

•All  numerical  results  shown  in  these  and  other  figures  are  for  an  ideal  gas  with  >  1.4. 
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second -order  potential  theory  due  to  Tobak  and  Wehrend  (Reference  13);  2)  the  unsteady 
flow  theory  due  to  Zartarian,  Hsu  and  Ashley  (Reference  14);  and,  3)  the  Newtonian 
impact  theory.  In  order  to  obtain  the  form  of  results  shown  in  Figures  29  and  30  from 
the  results  of  Zartarian,  et.  al. ,  it  was  necessary  to  substitute  expansions  of  the  form 


o  (t  -X/uj  *  a(t)  -  a  <t)~+  .... 


in  Equation  (38)  of  Reference  14. 


Figure  29  shows  a  fair  agreement  between  the  present  results  and  potential  theory 
for  a  10-degree  cone  at  low  Mach  numbers.  The  agreement  for  a  20-degree  cone 
(Figure  30)  is  not  so  good,  however.  At  higher  Mach  numbers,  the  impact  theory 
does  not  predict  Cn^q  and  as  well  as  might  be  expected  based  on  results  for 
CNa  •  The  agreement  with  the  theory  of  Reference  14  is  good  over  limited  Mach  No. 
ranges  for  both  the  10-degree  and  20-degree  cones  but  the  two  sets  of  results  tend  to 
diverge  at  both  high  and  low  Mach  numbers. 


Figure  31  shows  a  comparison  of  numerical  results  with  experimental  data  for  a 
10-degree  cone  at  Mach  10  over  a  range  of  center  of  gravity  locations.  The  data  also 
cover  a  range  of  Reynolds  numbers.  The  agreement  between  theory  and  experiment  is 
quite  good  except  for  the  highest  Reynolds  number.  Figure  32  shows  the  experimentally 
determined  variation  of  the  dynamic  stability  parameter,  Cm  +  Cm.  ,  with  Reynolds 
number  and  gives  some  indication  that  tne  disagreement  just  noted  may  be  due  to 
boundary  layer  transition. 
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Figure  33  shows  a  comparison  of  theoretical  results  and  experimental  data  for  j 
20-degree  cone  at  Mach  8.  The  data  are  for  varying  angle  of  attack  but  show  very  little 
influence  of  angle  of  attack.  The  agreement  between  theory  and  experiment  at  zero 
angle  of  attack  is  seen  to  be  quite  good. 

Figure  34  shows  the  theoretically  predicted  effect  of  frequency  of  oscillation  on  the 

dynamic  stability  parameter,  Cmn  +  Cm*  >  for  a  10-degree  cone  in  a  forced -oscillation 

q  ot 

wind  tunnel  experiment.  The  values  of  the  constants,  K  ,  K 2>  etc. ,  indicated  on  the 
Figure  34  are  related  to  the  unsteady  force  coefficients,  ~  ,  etc.  The  result 

is  that  there  is  no  detectable  effect  of  frequency  over  the  range  of  frequencies  en¬ 
countered  in  wind  tunnel  testing.  This  is  not  in  agreement  with  experimental  observa¬ 
tions  which  show  a  fairly  large  effect  of  frequency  (Reference  15).  The  effects  of 
frequency  appear  to  be  a  strong  function  of  Reynolds  number. 
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Figure  3.  Perturbations  In  the  Velocity  Component,  u 
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Figure  6.  Perturbations  in  t&e  Velocity  Component,  v 


u 


Figure  8.  Perturbations  In  the  Velocity  Component  ,  w 


ot|DS 


10.0  10.5  U.O  11.5  12.0  12.5 

*  ~  DEGREES 


Figure  12.  Perturbations  in  the  Entropy 
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10.0  10.5  11.0  11.5  12.0  12.5 

w  -  DEGREES 

Figure  14.  Perturbations  in  the  Stagnation  Enthalpy 
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Normal  Force  Derivatives 
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Figure  17.  Dynamic 


Dynamic  Normal  Force  Derivatives 


e  21.  Dynamic  Normal  Force  Derivatives, 


0  0.2  0.4  0.6  0.8  1.0 

x 

_££ 

L 


Figure  22.  Dynamic  Moment  Derivatives,  and 
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Figure  30.  Comparlaoa  of  Theoretical  Remits 
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Figure  31.  Comparison  of  Theory  and  Experiment 
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Figure  34.  Theoretical  Evaluation  of  Frequency  Effects 
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APPENDIX  n 


EQUATIONS  FOR  THE  FIELD  AT  ZERO  YAW 
The  following  equations  determine  the  unperturbed  field  at  sero  yaw: 
the  continuity  Equation, 

d  Po  vo  sin  ^ 

- - - +  2  p  .  u  8 in  to  “0; 

du>  o  o 

the  condition  of  irrotationality, 

d  Uq 

db!  ~  V0 

the  state  Equation, 

P0  =  P  V>o  So>i 

the  integrated  energy  Equation, 

S_  =  constant 
o 

and,  the  Bernoulli  Equation 

||2  ^  y2 

h  (P.») +  — g —  =  const- 
The  first  of  these  equations  can  be  re-written: 


and  la  uaed  in  deriving  Equation*  (5. 1),  (5.4),  and  (5. 6).  These  equations  are  used  to 
replace  derivatives  of  the  flow  field  quantities  at  aero  yaw  by  equivalent  expressions 
in  terms  of  the  flow  field  variables  themselves. 


APPENDIX  m 


EXPRESSIONS  FOR  THE  NORMAL  FORCE  AND  MOMENT 


Integration  of  the  perturbation  description  of  the  pressures  on  the  surface  of  the  cone 
gives  the  following  expressions  for  die  normal  force  (negative  of  the  force  In  the  y  direc¬ 
tion)  and  the  moment  In  the  z  direction  about  the  center  oS  gravity. 
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*  E 
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x  N, 
eg  l,n 

L  tan  x. 
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where: 


N 


_  PM 


l,n 


COB  UU 
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.  m  .  n+1 
(n+2)  sin  ou 
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2,  n 


2,n 


^ps  UU 
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n+2 

(n+3)  sin  iju 
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The  normal  force  derivatives  defined  in  the  symbols  section  are  related  to  the 


N,  „  as  follows: 
J.n 
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